Oscillatory settling in wormlike-micelle solutions: bursts and a long 
time scale 
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We study the dynamics of a spherical steel ball falUng 
freely through a solution of entangled wormlike micelles. 
If the sphere diameter is larger than a threshold value, 
the settUng velocity shows repeated short oscillatory bursts 
separated by long periods of relative quiescence. We 
propose a model incorporating the interplay of settling- 
induced flow, viscoelastic stress and, as in [PRE 66, 
025202(R) (2002) and PRE 73, 041508 (2006)], a slow 
structural variable for which our experiments offer inde- 
pendent evidence. 

'A sphere falling through a Newtonian fluid at low Reynolds 
jiumber Re^ reaches a steady terminal velocity through the in- 
terplay of gravity, buoyancy and viscous drag. This simple 
. .system, and its extensions to unsteady high Re flow, are well 
"understood and form the basis of falling ball viscometers^. 
In non-Newtonian fluids unsteady settling is observed ^-^ even 
at low Re. Mollinger et al. ^ observed periodic stick-slip set- 
.tling in the beginning of the fall and investigated the effect of 
wall confinement on the baseline settling velocity and the fre- 
.quency of unsteady events. Jayaraman et alA found a transi- 
tion from steady to irregular unsteady settling, remarked that 
the relevant control parameter was the shear rate, y = vi,/d 
(vi, & d being the baseline velocity and sphere diameter re- 
.spectively), and emphasized the importance of the negative 
wake behind the particle^'^. It was suggested'* that the oscil- 
lations were a result of the formation and breakup of flow- 
induced structures, known to arise in concentrated micellar 
.systems ^"^. 

In this paper we investigate the sedimentation of spherical 
metallic particles through an entangled wormlike-micelle so- 
lution^ of the surfactant CTAT (Cetyl Trimethyl Ammomium 
p-Toluenesulphonate) in brine. We report here for the first 
time that spheres larger than a critical size not only undergo 
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unsteady motion, but also show sustained, repeated bursts 
of oscillations superposed on a constant baseline velocity as 
shown in Fig. [T] The intra-burst oscillation period Zfast ^ 
d/vi,; the time between bursts isiow is at least an order of mag- 
nitude larger We account for these results through a model 
incorporating shear-induced coupled oscillations of viscoelas- 
tic stress components and a slow structural relaxation^SJi giv- 
ing rise to isiow We suggest that equilibrium fluctuations, i.e., 
without imposed flow, should show signatures of the structural 
variable. 

Experiments were conducted on two different samples. Gel 
A comprising 2 wt% CTAT + 100 mM NaCl and Gel B wifli 
2.2 wt% CTAT + 82 mM NaCl. Rheological measurements 
were carried out on a Paar Physica MCR 300 Rheometer in the 
cone-plate geometry. The frequency sweep and flow curve of 
Gel A give an extrapolated zero-shear-rate viscosity of 17.8 Pa 
s at 25°C, onset of shear thinning at a shear rate Yc ^ 0.4s^'. 
The linear viscoelasticity is Maxwellian with a relaxation time 
CO^/ ~ 1.34 s. The corresponding numbers for Gel B are 
8.5 Pa s, 1.5 s^' and 0.97 s respectively. These systems are 
in the parameter range where shear banding instabilities and 
spatiotemporal rheochaos are expected above a critical shear 



rati 



, 12-15 



A transparent cylindrical tube of diameter D = 5.4 cm and 
height 40 cm was filled with the wormlike-micelle solution. 
Steel balls of diameters d from 3 mm to 8 mm were introduced 
one at a time into the solution with zero initial velocity using 
an electromagnet positioned accurately to ensure that the balls 
sediment along the axis of the tube. Images of the falling balls 
were recorded using a Citius CI 00 Centurio camera capable 
of capturing 423 frames per second at its full resolution of 
1280x 1024 pixels. Image analysis was done on the captured 
frames to extract information about the position and veloc- 
ity of the ball as a function of time. Experimental runs were 
separated by at least 30 minutes in order to ensure complete 
relaxation and homogenization of the fluid prior to each run. 

We present first our results for Gel A. Fig. [Tlshows velocity- 
time plots of various diameters d. For d = Q.Acm (Fig. [IJa)) 
the ball quickly achieves a steady terminal velocity as ex- 
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Fig. 1 Velocity as a function of time: For Gel A; (a) d=0.4 cm, (b) 
d=0.56 cm and (c) d=0.7 cm. In (d), the velocity variation marked in 
(c) has been expanded. For Gel B: (e) d=0.5 cm, (f) d=0.7 cm and 
(g) d=0.8 cm. In (h), the velocity variation marked in (g) has been 
expanded. 



pected in the Stokesian limit for a Newtonian fluid. For 
d = 0.56cm (Fig. [Ub)) a strong oscillatory transient as in 
noted by King and Waters ^^ is seen before terminal velocity 
sets in. For d = 0.70cm (Fig. [Itc)), oscillations are seen to 
occur in repeated bursts on top of a constant baseline velocity 
Vh, and are sustained throughout the fall of the particle. The 
time between bursts is 4 to 7 sec, at least an order of magni- 
tude larger than the oscillation period ^^ 0.3 s within the burst. 
Fig- Ed) is an enlarged version of the bursts highlighted in 
Fig. ]Ilc). Within each burst the oscillations have an expo- 
nentially damped sinusoidal form; the total number of such 
bursts before the ball reaches the bottom is too small for us to 
conclude whether their timing is truly periodic or merely dis- 
tributed naiTowly about a well-defined mean. The figures also 
show that the oscillations are occasionally large enough to re- 
sult in a small negative velocity; i.e., the ball sometimes jumps 
up against gravity. Visual observation shows a strong negative 
wake behind the ball during the burst. Gel B shows similar 
behaviour as seen in Figs. [Ue-h). Interestingly, it should be 
noted that the baseline velocities scale as d^, consistent with 
Stokes's Law. Jayaraman et al.'^ remarked that oscillations set 
in at velocity gradients v^/d close to the frequency where G" 
(loss modulus) begins to increase. We do not observe such a 
correspondence. 

Fig. 2 shows the effect of wall confinement for a ball of 
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Fig. 2 The motion of 0.7 cm ball in CTAT (2wt%)+100 mM NaCl 
dropped in tubes of various diameters. Note the increase in the 
baseline velocity, Vj, in thinner tubes (diameter D) followed by 
decreasing time-period {T) of burst oscillations. 



diameter d=0.7 cm. A comparison of baseline velocities Vf, in 
tubes of diameter D = 5.4 cm, 3.2 cm and 2.0 cm shows that 
Vj, increases as D decreases, an effect opposite to that known 
in Newtonian fluids^. The effect is observed for smaller di- 
ameters as well. We find that within a burst the time-period 
T (as marked in Fig. |2]i decreases with decreasing D. Local 
thinning caused by flow alignment of micellar structure along 
the walls, expected to be more pronounced in narrower tubes, 
is a likely explanation for these observations. 

In order to account for the observation of bursts and oscil- 
lations, we now offer a simple theoretical model based on the 
interplay of the orientation field of the micelles and flow gen- 
erated by the settling particle. We emphasize that our intent 
is to offer a rational account of the observed phenomena, not 
to reproduce them in full detail. Further, we will show that 
the model, although schematic in nature, nonetheless leads 
to other predictions that can be tested in independent exper- 
iments. We start by summarizing some features of the be- 
haviour shown in Fig. [TJc). The phenomenon of primary in- 
terest is the observation of rapidly oscillatory bursts repeated 
at fairly regular long intervals. Each burst is damped in a 
roughly exponential manner, with a time constant noticeably 
smaller than the interval between bursts, i.e., each successive 
burst emerges from a flow that has become nearly quiescent. 
The baseline settling speed vq ~ 1.15 cm/s, with oscillations 
up to 12 cm/s. The resulting shear-rate at the scale of a par- 
ticle (diameter d ~ 0.7 cm) is y = v^/d 2± 1.64 to 17.1 s^\ 
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corresponding to a timescale y^' ~ .06 s to 0.6 s for the ac- 
cumulation of a strain of unity. The micellar solution has a 
measured zero-shear-rate viscosity of r|o — 18.0 Pa s, thinning 
under shear to about 1.4 Pa s, and a viscoelastic or orienta- 
tional relaxation time x,, ~ 1 .34 s. The Reynolds number at the 
particle scale is thus at most 6 x 10^^. The settling particles 
are steel (density p ~ 8 g/cm^) and hence would reach termi- 
nal velocity, if the viscosity were constant, in a time of order 
T, = p(i^/18r|o — 1.2 to 15.6 millisec depending on whether 
one takes the zero-shear or thinned value for the viscosity. The 
observed period within the bursts is Ti ~ 0.36 s, and the pe- 
riod between bursts is 72 ~ 9 s. Ti and Ty are comparable 
to Y^\ T2 ^ 71, Ty, and T,- is much smaller than all of these. 
Both acceleration and advection can thus be neglected in the 
Navier-Stokes equation, and we can work in the limit of Stoke- 
sian hydrodynamics, keeping track of stresses arising from the 
complex nature of the micellar solution. We now present our 
model for the coupled dynamics of the particle and the micel- 
lar medium. 

Consider a particle settling under gravity in a wormlike- 
micelle solution with velocity field u(r) at location r, with 
solvent viscosity r| and a contribution o'™' to the stress ten- 
sor from the conformations of the micelles. The steady Stokes 
equation 

-VP + W-V-o^"'^ (1) 



rjV^u 



balances viscous, pressure, and conformational stresses 
against the gravitational force from the buoyant weight W of 
the particle. Inverting dU and choosing the particle's centre of 
mass as origin, the settling velocity 



(2) 



v = vo+ d'rH{r)-V-a^"'>{r) 



where H is the incompressible Stokesian hydrodynamic prop- 
agator and vo is the Stokes settling velocity in the absence 
of contributions from micelle conformations, and the second 
term on the right in ([2) contains modifications of the settling 
speed due to micellar stresses around the particle. The dynam- 
ics of o''"' 



Da^'"'> , ^ 1 

-^ = XoA + Xi[A.o("V--" 

Dt T 



(ni) 



(3) 



combines flow-orientation coupling and relaxation, familiar 
from a variety of contexts including general rheological mod- 
els, rheochaos, and nematic hydrodynamics ' UZ"i2. Here 
D/Dt is a time derivative comoving and corotating with the 
fluid, A is the symmetric deformation-rate tensor, and the sub- 
script ST denotes symmetric traceless part. 

It is useful to work with a simple model extracted from 
Q and ([3]). From ^ we suggest that a component p of 
the micellar stress, presumably originating in micelle orien- 
tation around the particle alters the settling speed: v — V() + p 



where vq is the speed in the absence of contributions from the 
micellar stress. The form of (|3) suggests the schematic dy- 
namics p ^ v + vp + vq — p + ..., where production by shear 
competes with relaxation to local equilibrium. As in the 
Johnson-Segalman model '^ or nematic hydrodynamics'^— 
shear, through v, naturally couples p to other components 
q, which in turn have a similar dynamics. The ellipsis in- 
dicates terms of higher order in p and q, and we have as- 
sumed that the local shear-rate is proportional to the settling 
speed V. The structure of these equations raises the possibil- 
ity of oscillations with a timescale set by the shear-rate. Such 
shear-induced mixing of micellar stress components is very 
likely the mechanism for the oscillations within a burst, as 
their frequency, as seen from the numbers reported above, lies 
right in the middle of the range of particle-scale shear-rates y. 
The timescale associated with the repeated appearance of the 
bursts is far longer, as remarked above, and must arise from a 
distinct process. We propose that it enters through the mem- 
ory mechanism introduced by Gates and coworkers '"•" that 
we incorporate and explain in detail below. Assembling these 
ingredients we replace the dynamics of the micellar stress ten- 
sor ([Sj by the effective model 



p ^Xqv + Xi vp + X2vq — Z,p 



sw 



np,q), 



(4) 



q^fjQV + fJivp+fJ2vq-yq-(lip,q), (5) 

and (|2|i by the schematic relation v = VQ + p already mentioned 
above. The terms 2' and Q, in (IH, ^ incorporate nonlinear- 
ities that would arise naturally in a local thermodynamic ap- 
proac h'^'^" where p and q are related to the orientational or- 
der parameter. The coefficients X,-, /j,- are flow-orientation cou- 
plings, corresponding to similar terms in (|3) and best under- 
stood from nematic hydrodynamics^— . We have allowed 
the relaxation of p to occur through two channels: locally in 
time via the term Cp in (|4]i, and through the intervention of a 
structural variableRiSJi^ 



wit) 



dt e 



'J'-t'W 



Pit') 



determined by the past history of p. Equivalently, 

w = T{p — w), 



(6) 



(7) 



whence it is clear that w becomes identical to p only for 
r — ;> 00. We see formally that w and p stand in the same 
relation as position and momentum for an oscillator, even 
in regard to the relative signs of the coefficients of w in (|4|i 
and p in (jT). We emphasize that this relation has nothing to 
do with real inertia: our equations are meant to describe the 
viscosity-dominated regime. Equation d?) says, plausibly, that 

* In principle q could also relax through such a memory mechanism. For sim- 
plicity we have not explored this possibility. 
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Fig. 3 Oscillations in the settling speed from (|4j, (|5}, Q for several 
a as indicated. Note the two distinct timescales, within and between 
bursts. 



an imposed nonzero local alignment p shifts the equilibrium 
value of an underlying structural quantity such as the Maxwell 
viscoelastic time of the material^SJi or the micelle length^^. 
As in those works, we assume that the parameters F and s 
of this mechanochemical coupling are properties of the un- 
sheared system; our observations are consistent with F <C y, 
the particle-scale shear rate. 

The parameter space of (|4]l, (|5]l and dTjl is vast, and our lim- 
ited exploration finds regimes of spontaneous oscillation. Fig. 
[3] shows p (which is v apart from a constant) as a function 
of time for the model equations p—A\ —A^p —A^^q ~A/s,w + 
A^p''- -A(,p^, q^a(Bi~B2p-q^B^p'^),w ^Y[p-^v% 
bringing in q gradually by stepping a ~ 0,0.1,0.15&1.0. 

Our model, constructed using physically plausible argu- 
ments, reproduces key features of the experimental observa- 
tion, namely, a rapid oscillation timescale clearly associated 
with the shear and a much slower oscillation of distinct origin. 
Further, if the contribution of q is diminished by reducing a 
only the slow oscillation of p against w is seen. Only when a 
is large are the rapid oscillations seen, and these disappear on 
a timescale much smaller than that associated with the slow 
oscillation. 

As noted in ^''^^ , equations of the general form (|4]i - (|2) are 
widespread in the neuronal modelling literature ^^'^^ . Among 
the most robust features of such models are oscillations with 
a bursty character, and the coexistence of two or more widely 
separated timescales. In the neurophysics literature the intro- 
duction of an additional slow "recovery variable" was moti- 
vated in part by detailed knowledge of the electrophysiology 
of the axon—. The success of our model in capturing essential 

twith parameter values Ai=3, A2=10, A3=l, A4=A5=6, A6=l, Bi=3, -82=12, 
B3=6 & r=0.006 



features of our experiment encourages us to suggest indepen- 
dent experimental tests for the coupled dynamics of micellar 
stress p and "structure" w in the absence of shear. If not for the 
structural variable w we should have a conventional Johnson- 
Segalman or nematogenic fluid at rest, whose dynamical re- 
sponse and spontaneous fluctuations are overdamped. How- 
ever, in (|4|i and (|7j the micellar stress p and the structural 
variable w act like the momentum and coordinate of an os- 
cillator. And precisely because we require w to be slow, the 
damping coefficient F ^ s, the coupled dynamics is in the 
underdamped regime, so their linear response to disturbance 
must be oscillatory. Further, if we augment the dynamical 
equations (HI, (|7]) with thermal noise, the resulting Fourier- 
transformed dynamic coiTelation functions of the system in 
the absence of the falling ball should show a peak at a fre- 
quency around that of the slow oscillation seen in the falling- 
ball experiments. This feature seems to be quite robust, relying 
only on the assumption that the parameters in the dynamical 
equation for the structural variable w depend negligibly on the 
shear rate. This property is shared by the models presented 
by Gates et al.^^^ and emerges naturally from microscopic 
kinetic models^. However, it is conceivable that the kinetics 
in our system is such that w becomes slow only under strong 
imposed shear, in which case the equilibrium, i.e., unsheared, 
fluctuations of the coupled p — w dynamics would be over- 
damped. A search for such oscillatory response or correla- 
tion in equilibrium experiments would help resolve this issue. 
This in turn requires a dynamical probe of local anisotropy; 
a possible scheme might be to suspend droplets of a strongly 
flow-birefringent fluid in the wormlike-micelle solution, and 
measure the power spectral density of the spontaneous fluctu- 
ations in the light intensity seen through crossed polars. 

We emphasize that our simple model aims primarily to cap- 
ture the bursts and the two disparate time-scales in the prob- 
lem. Within the limitations of this approach, we have sug- 
gested possible independent tests of the model. Further ex- 
periments with different viscoelastic parameters and particle 
shapes should lead us to a more quantitative theory, at the level 
of reduced models like (HJi - (|7]l or a more complete hydrody- 
namic treatment. A better understanding, perhaps in terms of 
micelle kinetics^, of the slow structural variable in our system 
is of particular importance. 

To summarise: we have isolated the mechanism for the ob- 
served oscillatory descent of a sphere through a wormlike- 
micelle solution. The oscillations produced by (|4]i-(|7]i show 
two broadly different timescales, as in the experiments, al- 
though the detailed form of the oscillations is not identical to 
those observed. One timescale is associated with the shear, 
the other with the coupled dynamics of micellar stress and 
a structural parameter such as micellar length or equilibrium 
Maxwell time. We argue that this second timescale is likely to 
appear in the unsheared system in independent measurements 
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of the autocorrelation of spontaneous fluctuations of stress or 
orientation, and suggest experiments to test this possibiHty. 

Flow-birefringence studies ^^'^^ focusing on the velocity 
and orientation fields of the medium around the falling ball 
will offer further insight into the detailed mechanism underly- 
ing the phenomena presented in this work. We note that the 
birefringence contrast as seen by Rehage et al. ^^ in studies on 
solutions similar to ours, is small. Such measurements on our 
system will therefore require a careful high resolution study, 
which we defer to later work. 

Numerical studies of the complete hydrodynamic equations 
for orientable fluids, studying the effect of the gravitational 
settling of a sphere through the medium, should see these cu- 
rious oscillations, and will allow us to converge on the correct 
effective model of the general type of (|4]i - (|7]i- 
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